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We study electromagnetic field propagation through a planar three-waveguide coupler 
with linear gain and loss, in a configuration that is the optical analog of a quantum 
VF -symmetric system, and provide its closed-form analytic propagator. At an spe¬ 
cific propagation length, we show that the device provides all-optical amplitude to 
phase modulation with a n modulation range, if an extra binary phase is allowed in 
the reference signal, as well as phase to amplitude modulation, with an amplitude 
modulation range that depends linearly on the gain-to-coupling ratio of the system. 
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Optical communications are a necessity of modern life . 12 All-optical systems based on 
monolithic photonic integrated circuits may prove a compact, robust and ultrafast alter¬ 
native with low-dissipation to some optoelectronic devices )- 3 For example, since the report 
of the optical directional couplers with attenuated waveguides 4 and the nonlinear coherent 
coupler,® three-waveguide nonlinear directional couplers, based on the Kerr nonlinearity, 
have been shown to provide all-optical spatial switching 6 12 and logic gates.®^® 

Here, we are interested in the planar three-waveguide coupler with a gain mechanism 
different from the Kerr nonlinearity. In particular, we are going to study waveguides with 
complex effective refractive index, in other words, complex effective first-order susceptibili¬ 
ties, such that we have effective linear gain or loss in each waveguide. Such a system may 
be described by coupled-mode theory 4 16 18 in the form of the differential equation set, 



= n 0 S 0 (z) + g£i(z), 

( 1 ) 


= ni£i(z) + g [£ 0 (z) + S 2 {z )\, 

( 2 ) 


= n 2 S 2 (z) + gSi(z ), 

( 3 ) 


where the complex held amplitude at the j-tli waveguide is given by £j(z), the constant 
effective refractive index is rij = ujj + Efj , with real parameters ujj and 7 j, and the constant 
effective coupling between waveguide modes is g. While it is possible to deal with this 
coupler, we will focus on the so-called VT -symmetric optical systems. 19 21 In our device, 
that translates into effective refractive indices of the extremal waveguides that are complex 
conjugates of each other and that of the central waveguide is real. This model can be realized 
by engineering a three-waveguide coupler where the real part of the effective refractive 
index is ujq = uj 2 = uj and the gain-loss at each waveguide fulfills 71 = (70 + 72) /2 and 
7 = (70 — 72) /2, such that a solution of the form, 


£j{z) = e^+^Ejiz), 


( 4 ) 


delivers an effective coupled-mode differential set, 

d 


- i^E 0 (z) = (u - ivi + i'y) E 0 (z) + gE x (z), 
dz 

-* T~ El(yZ ^ = 9 1 E ° + e<2 ( z ^ > 

-i~T E 2 {z) = (w - uji - i'y) E 2 (z) + gE x {z). 
dz 
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( 5 ) 

( 6 ) 
( 7 ) 




Amplitude to phase modulation is another key element in communications and, to the 
best of our knowledge, an all-optical solution involving TT-symmetric three-waveguide cou¬ 
plers has not been reported in the literature. For that reason, in the following, we will discuss 
propagation through the three-waveguide PT-symmetric coupler in the case of waveguides 
with identical real part of the effective refractive index, provide a closed-form analytic prop¬ 
agator for the device and explore its use as an amplitude to phase and phase to amplitude 
modulator. 

For the sake of simplicity, we will consider waveguides with identical real part of the 
effective refractive index, uj± = lu, such that we obtain the coupled-mode set, 

-i^E(0 = ME(0, ( 8 ) 

where we have defined a scaled propagation distance, £ = A z, in order to have a single control 
parameter in the form of the ratio between the imaginary part of the effective refractive index 
and the coupling parameter, £ = 7/A, henceforth named gain-to-coupling ratio, and used 
the shorthand notation, 


H = 


if 1 0 \ 
10 1 
0 1 -if y 


e{ 0 = 


( E 0 {Q\ 

m), 

{ £2(0 J 


( 9 ) 


for the mode coupling matrix and the complex field amplitudes vector, in that order. It is 
straightforward to show that this mode coupling matrix has two nonzero eigenvalues, ±h2, 
with 


ft = v / 2 r e (10) 

The eigenvalues will be real when the gain to coupling ratio fulfills £ < a/2, completely 
degenerate for £ = y/ 2 , and complex if £ > \/ 2 . In other words, VT -symmetry is broken in 
the latter case. 

The mode coupling matrix, H, does not depend on the scaled propagation distance, £, 
thus, we can calculate the propagated complex field amplitudes asp^ 

E( C) = U(C)T( 0 ), U(C) = e®/ (11) 


where the input complex field amplitudes are collected in the vector 22 ( 0 ). The propagator, 

U(£), can be calculated from its Taylor series, 

00 (iCV 

u(0 = ( 12 ) 

j =0 T 
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( 13 ) 


f 1 + ^sinnCH + i (cos^C - 1) H 2 , £ ^ V2, 

\ 1 + *CH- K 2 H 2 , e = V^’ 

where the identity matrix is given by 1 and we used the Cayley-Hamilton theorem,® a 
square matrix satishes its own characteristic equation, in this particular case, 

H 3 = H 2 H. (14) 


After some algebra, it is possible to provide the propagation as a response to an impulse 
function, that is, the complex held amplitude at the n-th waveguide given that the input 
complex held amplitude impinged at just the m-th waveguide, 


4 m) (0 = U„, m (C), n,m = 0,1,2, 

v m+n 


a 


cr 



2\ 

M r 



2 ih 

i 


\ n ) 


—771, —77, —2, — 


ot 0 
a 2 \ 


with Gauss hypergeometric function, 2 i r i(a, b, c, d),® a nd auxiliary functions, 


= f &//3 2 , £ ^ y/2, 

*° \2/(V2-C) 2 ,{ = V2, 

I iy/2/fi, £ / y/2, 

a = < 

[iC/(V2~ C), £= A 

/3 = ficos -f2£ — £sin -f2£. 


(15) 

(16) 


(17) 

(18) 
(19) 


Thus, even with amplihcation and losses, we expect periodic oscillations of the complex held 
amplitudes in the VT -symmetric regime, £ < y/2 with a period equal to 27rfW 1 , Fig. [lj that 
will disappear as the VT -symmetry is broken, £ > y/2, and gain dominates, Fig. [l} We 
conducted both a numerical solution of the differential equation set and propagation with 
the analytic results and obtained good agreement between both methods. 

In the following, we will focus on just the VT -symmetric regime, 0 < £ < y/2, because 
the device always behaves like a directional amplifier in the broken symmetry phase, Fig. [2j 
We want to bring forward that the response to impulse function, Fig. [lj reveals that, at a 
certain scaled distance, 


C/ 


arccos (£ 2 — 1) + 2m r 

n 


77 = 0,1,2, ..., 


( 20 ) 
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FIG. 1. (Color online) Squared response to impulse function, (Fn"' 1 1 2 , at the zeroth (dotted blue), 
n = 0, first (solid black), n = 1, and second (dashed red), n = 2, waveguides for an initial held 
impinging at the (a) zeroth, m = 0, (b) first, m = 1, and (c) second, m = 2, waveguide of a coupler 
with gain-to-coupling ratio £ = 0.5 that keeps PT-symmetry. 



FIG. 2. (Color online) Same as Fig. [ 2 ] but for a gain-to-coupling ratio £ = 1.25\/2 that breaks 
VT -symmetry. 
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the output complex field amplitudes are quite simple, 


£o(C/) = -M o), (21) 

m f ) = — -^i(o)+*2^2(0), (22) 

MCf) = -Eq(0) + 2£ [^(0) + £E 2 {0)] . (23) 

If we consider initial complex field amplitudes, 

E 0 { 0) = Aoe** 1 , Ex(0) = Ae #1 , £ 2 (0) = 0, (24) 


with real amplitudes and phases, A 3 > 0 and 0 < <fi 3 < 2n, the fields at the output are given 

by, 


Eo(Cf) ~ 0, 

E^Cf) = i ie i( ^ + b, 

E 2 (C f ) = Ae^ 0+7r) + 2^A 1 e i ( rt>1+ ^. 


(25) 

(26) 
(27) 


Thus, it is possible to convert amplitude and phase differences between the signal, £o(0), 
and reference, E'i(O), input fields to phase and amplitude differences at the ontpnt field, 


\E 2 (Cf) \ = \[j§ + ^ 2 A\ - ^AqAi sin 5, 
A 0 sin 0 O - 2£Ai cos fa 

arg [E 2 {Q)\ = arctan 


(28) 
(29) 

Aq cos 0o + 2 £.Ai sin 0 i 

with the phase difference between signal and reference given by 8 = 0 O — <Pi- We will take 
the argument function as modulo 2n to have a phase in the range [0, 2n). Note, we will 
assign a null phase value to zero fields. 

Let us consider first the case of input fields of proportional amplitudes, 

E 0 (0) = yd -Ai. E x ( 0 ) = ±Ai, E 2 ( 0 ) = 0 . (30) 

Note that we have added a binary phase on the reference field, £i(0), that will help ns in 
the future. The output field amplitude and phase are given by the following, 

|S(C/)I = (31) 

arg [E 2 (Q)} = arctan t2£t?, (32) 


where the phase is a function of the reference-to-signal amplitude ratio, 

Ai 


V = 


i-Aj 


(33) 
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FIG. 3. (Color online) (a) Squared field amplitude, \Ej(Q)\ 2 , and (b-c) phase, arg Ej(Q), with 
j = 0 (dotted blue), j = 1 (solid black), and j = 2 (dashed red), as a function of the reference field 
(a) E'i(O) = A\ and (b) E\ (0) = —A\ for a PT-symmetric amplitude to phase modulator with 
gain-to-coupling parameter £ = 0.5. 

For any given gain-to-coupling ratio, each of the two binary phase options provides a 7 t/ 2 
modulation range but in different phase regions, 

= [arg£ 2 (C/)] max - [arg£ 2 (C/)] min , 

7T-7T/2, (0) — Ai i 

37t/2 — 7r, ^i(O) = —Ai, 

7T 

= 2 ’ 

Thus, the introduction of the binary phase allows for continuous modulation of phase in a 7r 
range. Figure [3] shows an example of amplitude to phase modulation with a gain-to-coupling 
ratio £ = 0.5 that provides an stable output field amplitude, Fig. [3](a), for all the phase 
modulation range, Fig. [3](b) and Fig. |3](c). 

The converse situation is provided when the input fields have the same amplitude but 
different phase, 



(34) 

(35) 

(36) 


So(0) = e“A u £,(0) = A u E 2 ( 0) = 0. 


(37) 
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FIG. 4. (Color online) (a) Squared field amplitude, \Ej(Q)\ 2 , and (b) phase, arg Ej{C,f), with 
j = 0 (dotted blue), j = 1 (solid black), and j = 2 (dashed red), as a function of the initial 
phase difference, 5, between signal, Eq( 0) = e tS / \/2, and reference, E\(0) = l/\/2, fields for a 
FT-symmetric phase to amplitude modulator with gain-to-coupling parameter £ = 0.5. 


Here the output field amplitude and phase are a function of the phase difference, 


\E 2 (C f )\ = ^l + 4£(£-l)sin5 A , (38) 

arg [T 2 (C/)] = arctan —— S ^. (39) 

cos o 

Thus, we are mapping the phase difference to both the amplitude and phase of the output 
field. The range of the phase to amplitude modulation can be measured by an equivalent of 
the interferometric visibility, 


I^VC/IILx - |£2(C/)IL. 

|£2(C/)IL* + |£ 2 (C,)IL„’ 

SfkAjvTi 

A$ + 4?Al' 


(40) 

(41) 


where we can see that the gain-to-coupling ratio plays a fundamental role. Figure [4] shows 
the amplitude and phase of the output fields as a function of the initial phase difference 
between signal and reference fields, for a gain-to-coupling ration £ = 0.5 with equal initial 
complex field amplitudes, A\ = A 2 = 1 / a/2- 

In summary, we have studied the propagation of electromagnetic fields through three 
coupled waveguides with gain. The real and imaginary parts of the couled mode effective 
refractive index of the waveguides are chosen to produce an optical analog of quantum VT - 
symmetry. We focus in the particular case of refractive indices with identical real parts and 
provide a closed form analytic propagator for the device. Then, we show that the device can 
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be used as an optical oscillator or directional amplifier inside and outside the PT-symmetric 
regime, in that order. 

Furthermore, we show that there exists a characteristic length that delivers a two-input, 
single-output, VT -symmetric device providing all-optical amplitude (phase) to phase (am¬ 
plitude) modulation. The amplitude to phase configuration provides a phase modulation 
range of 7t/2 that can be increased to 7r if a binary phase is introduced into the reference 
signal. The phase to amplitude mode provides an amplitude modulation range that depends 
linearly on the gain-to-coupling ratio of the waveguide array. 
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